Low Temperature Penetration depth of K-(ET)2Cu[N(CN)2]Br 

and /€-(ET)2Cu(NCS)2 
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We present high precision measurements of the penetration depth of single crystals of 
K-(ET)2Cu[N(CN)2]Br and k-(ET)2Cu(NCS)2 at temperatures down to 0.4 K. We find that, 
at low temperatures, the in-plane penetration depth (Ay) varies as a fractional power law. Ay ~ T'. 
Whilst this may be taken as evidence for novel pair excitation processes, we show that the data are 
also consistent with a quasilinear variation of the superfluid density, as is expected for a d-wave su- 
perconductor with impurities or a small residual gap. Our data for the interplane penetration depth 
show similar features and give a direct measurement of the absolute value, A_l (0) = 100 ± 20 pm. 



PACS numbers; 74.70.Kn, 74.25.Nf 

Compounds of the family k— (ET)2X have the highest 
transition temperatures of all organic superconductors 
known to date |l|. They have recently attracted consid- 
erable attention because of their similarity to the high 
Tc cuprates and the possibility that they may also have 
a non-conventional paring state B. The two materials 
studied here, K-(ET)2Cu[N(CN)2]Br (Tc - 11.6 K) and 
K-(ET)2Cu(NCS)2 {Tc - 9.6 K), are highly anisotropic, 
layered, extreme type II superconductors. As in the 
cuprates, the superconducting phase in these materi- 
als is in close proximity to an antiferromagnetic phase. 
Both antiferromagnetic spin fluctuations and a pseudo- 
gap have been detected in NMR measurements in the 
normal state ||]. Neither the underlying pairing mecha- 
nism nor the symmetry of the order parameter has been 
conclusively established. Although NMR [§,^, specific 
heat g and thermal conductivity measurements all 
suggest a non-conventional pairing state, results of pene- 
tration depth measurements have been inconsistent, with 
evidence for both conventional g,g] and non-conventional 
[[10|-p2[ behavior. However, none of these penetration 
depth measurements have been performed over a temper- 
ature range {T/Tc) and a precision, comparable to those 
in the cuprates p3| . In this Letter, we present mea- 
surements of both the in-plane Aii, and the interplane, 
A±, penetration depths in K-(ET)2Cu[N(CN)2]Br and 
K-(ET)2Cu(NCS)2 at temperatures down to 0.4 K. 

Our measurements were performed on single crystals of 
K-(ET)2Cu[N(CN)2]Br and k-(ET)2Cu(NCS)2 which 
were grown at Argonne National Laboratory. Details of 
the growth procedures have been given elsewhere El. 
Penetration depth measurements were performed using 
a 13 MHz tunnel diode oscillator |1^ mounted on a ■^He 
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and low 



refrigerator. The low noise level [-^ 
drift of the oscillator allows us to obtain high resolution 
data with a very small temperature spacing interval. The 
samples were attached, with a small amount of vacuum 
grease, to a sapphire rod which fitted inside the copper 
sense coil. The sense coil was calibrated using spheres of 
Aluminum. The sample temperature was measured with 
a calibrated Cernox thermometer attached to the other 
end of the sapphire rod. The samples were cooled slowly 
(0.1-1.0 K/min) to avoid introducing disorder Q. The 
entire cryostat was surrounded by a triple mumetal shield 
that reduced the stray dc field to less than 0.005 Oe. The 
RF probe field was estimated to be ~ 0.005 Oe. This very 
low measurement field is important as in these materials 
Hci is very low (especially for H parallel to the layers) 
and the vortices are very weakly pinned in the mixed 
state [0|8| . The complex vortex dynamics are a serious 
problem for extracting A(T) from mixed state magneti- 
zation measurements |19| , and have not been fully taken 
into account in some previous studies faJy]. 

For measurements of Aii , the RF field was directed nor- 
mal to the conducting planes thus ensuring that only in- 
plane supercurrents are excited CO]. We relate the mea- 
sured frequency shifts to changes in A by measuring the 
area of the faces parallel to the field and multiplying by a 
demagnetization factor, which was estimated by the in- 
scribed or exscribed ellipsoid approximation [pi| . The as- 
pect ratios of our crystals varied from approximately 1.0 
to 2.5 (the short dimension is normal to the planes), and 
so the demagnetizing correction 1/(1 — N) was between 
1.6 and 2.4. We estimate that this calibration is accurate 
to ~20%, which is in accord with the measured variations 
of AAii between samples. To determine Aj_ the field was 
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FIG. 1. Temperature dependence of the in-plane penetra- 
tion depth A|| of a crystal of K-(ET)2Cu[N(CN)2]Br. The 
inset shows the susceptibility over the full temperature range. 



directed along the planes. The strong anisotropy of these 
materials ensures that in this orientation the response is 
dominated by interlayer currents for our crystals' aspect 
ratios. 

In Fig. nl we show the temperature dependence 
of A|| for a crystal of K-(ET)2Cu[N(CN)2]Br. Our 
K— (ET)2Cu(NCS)2 samples show essentially the same T 
dependence (see later). The rapid variation of A||(T) at 
low temperature strongly suggests the existence of gap- 
less excitations. The similarities between these materi- 
als and the cuprates suggest that these excitations may 
originate from a d-wave pairing state. In a pure d-wave 
superconductor it is well established that A|| will vary 
linearly with T at the lowest temperatures. Although 
the data in Fig. |l| show no region where A|| (T) is strictly 
linear there are several factors which may intervene to 
mask this simple power law behavior. 

First, Kosztin and Leggett [^ have argued that at 
low temperatures, of order T^. ^||/A||, nonlocal effects can 
change the linear behavior to quadratic. Both present 
materials however, are in the clean, local limit, A|| 3> 
i S> ^11 |11[ and so this effect should only be important 
for T ^ 70 mK, which is well below our minimum tem- 
perature. Next, there is the important distinction be- 
tween AA(r), which is directly measured, and the su- 
perfluid density [p(r) = }?{^)I'>?[T)\ which can only be 
inferred with knowledge of A(0). In the d-wave model, 
the non-linear corrections to p.. (T) at higher temper- 
ature are considerably smaller than for A||(T). If p.. 
varies strictly linearly with T ^ p — 1 — aT /Tc, then 
A(r) = A(0)[1 + \a{T/T,) + ^a^{T/T,)^ + •••]. There is 
therefore invariably a quadratic component to A whose 
relative strength depends on a which in the d-wave model 
is inversely proportional to the angular slope of the en- 
ergy gap at the nodes, a = dAte)/deUd^ lH' ^ ^^^^'^ 
factor is the influence of impurities. Bonn et al. EJ] 







FIG. 2. In-plane superfluid density p,, = Aj](0)/A[| (T)^ cal- 
culated from the AA||(T) data in Fig. 1, for several values of 
A||(0). The thin lines are fits to the data with Eq. Q. The 
inset shows the same data over a wider temperature range. 



have shown that in YBa2Cu306.95 even small amounts 
of impurities can change the low temperature behavior 
of p{T) from linear to quadratic. Hirschfeld and Gold- 
enfeld 1 25 1 have explained how this behavior results from 
impurity scattering in a d-wave superconductor and have 
suggested the following formula to interpolate between 
the high temperature (pure) and the low temperature 
(impurity dominated) regimes 



Ap„(T) 






(1) 



here T* parameterizes the impurity scattering rate. For 
T* to be significant without Tc being significant de- 
pressed requires that the scatterers be unitary P5| , p4| . 
A residual small energy gap at the 'node', resulting from 
a mixed order parameter (e.g. d^^-y^ + i\f^s) results in 
similar behavior |2^ . 

In order to compare our data to the d-wave model we 
have plotted in Fig. g the calculated superfluid density 
pq . Muon spin relaxation measurements [ [ll]j2^ give 
values of A||(0) = 0.78 ± 0.2 pm and 0.69± 0.2pm for 
K-(ET)2Cu[N(CN)2]Br and k-(ET)2Cu(NCS)2 respec- 
tively. However, because of the very weak pinning and 
strong anisotropy of these compounds, extracting A(0) 
from /iSR measurements can be problematic p9| and 
therefore we show p calculated for A||(0) in the range 
0.8 to 1.8 pm. The thin lines in the figure show a fit to 
the data with Eq. H^. It can be seen that the fit is reason- 
able for all values of Ay (0) but is best for the values close 
to the high end of our range. For A(0) = 1.8pm there is a 
quasilinear region extending from roughly 1 K to 6 K, and 
the fit to Eq. | gives a = 1.2 and T* = 0.6 K. Lower val- 



ues of A|| (0) results in a lower value of T* but restricts the 
upper limit of the quasilinear region; for A||(0) = 0.8/zm, 
T* = 0.4 K and the curve breaks away from the linear fit 
at 1.7 K. These values of T* /Tc are roughly comparable 
to what is found for good quality YBa2Cu306.95 and 
other cuprates |ljj24,E3,p3| . If this low T turnover is in- 
terpreted as arising from a mixed order parameter then 
our T* values would imply that the minimum gap is ^ 
3% of the maximum ||2|] . 

For the lower values of A||(0) the shape of the P«{T) 
curve is unusual in that it has positive curvature. Al- 
though we know of no cuprates which show this behav- 
ior it has been predicted in some d-wave models [ pO[ . 
The shape of p^^iT) for A||(0) ~ 1.8/im is similar to 
that reported for the single layer cuprate, Tl2Ba2Cu06+5 
p9| . The values of a extracted from our fits are some- 
what larger than that found for YBa2Cu306.95 (a ~ 
0.6 [|l^,^) but are comparable to that found for 
TlaBaaCuOe+a Hi (a - 1.0). This would imply that 
if K— (ET)2Cu[N(CN)2]Br is a d-wave superconductor 
then it has a small value of dA{9)/d9\node- We conclude 
therefore, that our data are compatible with a d-wave or- 
der parameter once the influence of impurities (or a small 
residual gap) and a large value of a are taken into ac- 
count. Although the above analysis has been conducted 
for K— (ET)2Cu[N(CN)2]Br, essentially the same conclu- 
sions are reached for k-(ET)2Cu(NCS)2 as AA||(T) for 
the two compounds are very similar (see below). The 
major difference is that if the p., (T) curves are forced to 
lie on top of those in Fig. then the A||(0) values for 
K-(ET)2Cu(NCS)2 are ~35 % lower (the lower limit of 
A||(0) in Fig. || becomes ~0.52/im). 

We now consider an alternative interpretation of the 
data. Kosztin et al. [M| have recently proposed a BCS 
to Bose-Einstein (BCS-BEC) crossover theory for short 
coherence length superconductors in which finite momen- 
tum pair excitations in addition to the usual fermionic 
Bogoliubov quasiparticles can deplete the condensate. 
This theory has been used to explain the pseudogap be- 
havior in the underdoped cuprates, and gives rise to a T^ 
term in A below a characteristic temperature T^g. NMR 
measurements indicate that pseudogap features are also 
present in the ET compounds [Bl which leads us to con- 
sider this theory as a possible explanation for our data. 

In Fig. H we have plotted our AA||(r) data for two 
samples of k— (ET)2Cu[N(CN)2]Br and two samples of 
K— (ET)2Cu(NCS)2 versus T2 . It can been seen that the 
data follow this law very closely over almost a decade of 
temperature. A fit of the low temperature data to a vari- 
able power law, AAj|(r) = /ST", gives n = 1.50 ± 0.04, 
illustrating the robustness of the fractional exponent. 
Within the errors the slope (3 (not normalized to Tc) was 
the same for all samples {(3 — 0.06 ± O.Ol/xm K~3 ). The 
upper limit of the power law is ~ 4 K but is somewhat 
sample dependent. 
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FIG. 3. AA||(r) data for two sam- 

ples of K-(ET)2Cu[N(CN)2]Br (a, b) and two samples of 
K-(ET)2Cu(NCS)2 (c, d) plotted versus [T/Tc\^. The data 
have been offset for clarity. The inset shows the data for 
samples a and d on a log scale. 



In the BCS-BEC theory both the magnitude of the T 2 
term in A(r) and Tgg depend on material parameters. 
By modelling K-(ET)2Cu[N(CN)2]Br with a simple 
quasi-2D tight binding band structure and a moderately 
strong coupling constant, it can been shown pH l that 
A|| (T) follows a Ta power law up to '-^ Tc/3. The crucial 
assumption of this fit is that the order parameter should 
have s-wave symmetry. To get quantitative agreement 
requires that A||(0) ~ 3.2/im for K-(ET)2Cu[N(CN)2]Br. 
So although the theory explains well the power law dis- 
played in Fig. ^, quantitative agreement requires that 
A|[(0) for our crystals is around 4 times the /xSR esti- 
mates. 

Finally in Fig. we show data for the temperature de- 
pendence of the interplane penetration depth for a sam- 
ple of K-(ET)2Cu[N(CN)2]Br. As Ai(0) is large we are 
able to determine its absolute value directly from our 
measurements. We obtain Aj_(0) = 100 ± 20/im which 
is comparable to previous measurements ||l^,^,^. As 
shown in the figure, p±_ is also strongly temperature de- 
pendent at low temperature. The relative strength of 
the temperature dependence of A in the two directions 
is AAi/AA|| ~ 17 (between 0.4 K and 3 K). p^[T) like 
p.. (T) is close to being linear at low temperature. A vari- 
able power law [p_L (T) = 1 — fiT^] fits the data well and 
gives an exponent n = 1.2 ± 0.1 (the uncertainty reflects 
the sample and T fitting range dependence). A dirty 
d-wave fit [Eq. (1)] is slightly worse and gives values of 
T* ~ 0.5 K and a = 0.42 ± 0.07. 

In these strongly anisotropic materials it has been ar- 
gued that pj_ depends on both in-plane scattering pro- 
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FIG. 4. The inter- 

plane superfluid density [p±_ — \\{Q) / }?^{T)] for a crystal 
of K-(ET)2Cu[N(CN)2]Br. The inset is an expanded view of 
the low temperature data. 



cesses as well as interplane hopping amplitudes and so 
is not in general expected to have the same T depen- 
dence as P|| [gj|. In HgBa2Cu06+5 it was found that 
although A|| ~ T, A^ ^ T^ |3q|, whereas in clean sam- 
ples of Bi2Sr2CaCu208-f5 both A|| and Aj_ vary linearly 
with T |36|. The fact that our data show a close to lin- 



ear T dependence however, supports our conclusion that 
these materials have low lying excitations. 

In conclusion, our measurements give strong evi- 
dence for the existence of low lying excitations in both 
K-(ET)2Cu[N(CN)2]Br and k-(ET)2Cu(NCS)2. There 
are two competing explanations for our in-plane data. 
We have shown that the observed non-linear T depen- 
dence of All can arise from a d-wave (or other gapless) 
order parameter with a small value of dlS.[9) / d9\node^ and 
either a small residual gap (^ 3% of the maximum) or 
a small density of unitary impurity scatterers. An alter- 
native view is that these materials have a gapped order 
parameter and that the T dependence of Aii arises from 
incoherent pair excitations expected in short coherence 
length superconductors which exhibit a pseudogap [ pl[ . 

3 

Whereas the observed Ta dependence of Aj| arises natu- 
rally in this model, the literature values of A||(0) are sig- 
nificantly closer to those required for the d-wave interpre- 
tation. Resolution of this issue will require measurements 
of A[|(T) to lower T and more accurate determinations of 
A||(0). 
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